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Abstract 

In this paper, we present the exact solution to a one-dimensional, two- 
component, quantum many-body system in which like particles interact with 
a pair potential s(s + l)/sinh 2 (r), while unlike particles interact with a pair 
potential — s(s + l)/cosh 2 (r). We first give a proof of integrability, then derive 
the coupled equations determining the complete spectrum. All singularities 
occur in the ground state when there are equal numbers of the two compo- 
nents; we give explicit results for the ground state and low-lying states in this 
case. For s > 0, the system is an antiferromagnet/insulator, with excitations 
consisting of a pair-hole-pair continuum, a two-particle continuum with gap, 
and excitons with gaps. For — 1 < s < 0, the system has excitations consisting 
of a hole-particle continuum, and a two-spin wave continuum, both gap-less. 
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I. THE S-C MODEL 



We present the exact solution to a one-dimensional, two-component, quantum many- 
body system of considerable complexity. The two kinds of particles are distinguished by a 
quantum number a = ±1, which may be thought of as either spin or charge. The system is 
defined by the Hamiltonian 

1 d 2 

H = ~ 9^2 + E Vjk(Xj-X k ), (1) 

l<j<N UX j l<j<k<N 

where the pair potential is 

+ CTjCTfc 1 — OjO-fc 



v jk {x) = s(s + 1) 



(2) 



2sinh 2 (x) 2cosh 2 (x) 

We assume s > — 1. We call this the S-C model, for the sinh-cosh interaction. Thus for s > 0, 
like particles repel, while unlike particles attract. When like particles are near, the repulsive 
potential increases as 1/r 2 , while for large separations, both potentials decay exponentially 
with a decay length we take as our length scale, and hence unity. The potentials might 
usefully be thought of clS db screened 1/r 2 potential. 

This system was first introduced by Calogero who showed it to be integrable. Suther- 
land |2| soon afterward showed that the system could be exactly solved, and gave the solution 
for a single component system. He further showed the Toda lattice to be the low-density 
limit, and was able to take the classical limit to reproduce Toda's celebrated results, iden- 
tifying the particle-hole excitations of the quantum system with the soliton-phonon modes 
of the classical system. Our present solution for the two-component system exploits in a 
fundamental way the integrability of the system, so we first discuss this point. 



II. INTEGRABILITY 

For a classical system of N one-dimensional particles, Lax [Q, Moser |3j and Calogero || 
have shown that for certain potentials one can find two Hermitean N x N matrices L and 
A that obey the Lax equation dL/dt = i[A, L\. Thus L evolves by a unitary transformation 
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generated by A, and hence det[L — a;!] is a constant of motion. Expanding the determinant 
in powers of lu, we find N integrals of motion, i.e., det[L — uit] = J2o<j<N Jj(~~ LJ ) N ~^ 3 — 
1, . . . , N. Further, these integrals have been shown to be in involution, and thus the system 
is integrable. 

Under a reasonable assumption for the form of the Lax A and L matrices, Calogero 
has shown that the most general solution to the Lax equations is given in terms of the 
Jacobi elliptic function sn(x|m), and by a suitable choice of parameters, our Hamiltonian is 
included. 

Calogero has also demonstrated that if one replaces the classical dynamical variables 
with the corresponding quantum mechanical operators, det[L — cull] is well defined with no 
ordering ambiguity, and the quantum mechanical commutator \H, det[L — cull]| = 0. Thus, 
the Jj are still constants of motion. Finally, Calogero showed that [ det[L — cull], det[L — 
o/lL]J = 0, and thus the quantum system is also completely integrable. 

To sum up the situation: Our system is completely integrable. For the general classical 
system integrability tells us something concrete, namely that the motion in terms of action- 
angle variables is on a torus. However, for the general quantum system integrability seems to 
buy one almost nothing. The exception is for those special cases which support scattering; 
i.e., systems which fly apart when the walls of the box are removed. In these cases, in the 
distant past and future the Lax matrix L approaches a diagonal matrix with the momenta 
as diagonal elements, so that det[L — ut] = l\i<j<N(Pj ~ u )- Thus the individual momenta 
Pj are conserved in a collision, and hence, as emphasized by Sutherland [[|, the wavefunction 
is given asymptotically by Bethe's ansatz. Sutherland has exploited this fact to completely 
determine, in the thermodynamic limit, properties of systems interacting by such potentials 
including our system. We stress that no features of the proof of integrability are needed. 
All we need is to know it to be integrable, by whatever method. 

The proof of Calogero, however , is very difficult, and only briefly sketched in the litera- 
ture. For that reason, we now offer an alternative proof of integrability based on a method 
of Shastry ||. Let us write the Lax matrices as 
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A jk = S jk + (S jk - (3) 

Ljk = SjkPj + i(l - S jk )a jk , (4) 



where 

ajk = —s — - — — -coth(xj — Xfc) H — -tanh(xj — x k ) (5) 

Then if the two-body potential potential Vj k is given by Vjk = oij k + a'j k — s 2 , we find that 
the quantum Lax equation \H, L\ = [L,A] is satisfied. Here the first fancy commutator is 
a quantum mechanical commutator between operators, while the second commutator is an 
ordered matrix commutator, so that the equation above is really iV 2 equations of the form 
[H, Ljk] = Y^KKNiLjiAik — AjiLik). This potential, however, is exactly the potential for our 
system. 

We now observe that the Lax A matrix has the following very important property. 
Defining a vector ( with Q = 1, we see A( = (^A = 0. This allows us to construct 
constants of motion by I n = (^"C since 

J B |= Clff.LniC (6) 

= C f E {L j \H,L\L N -^}C (7) 

0<j<N-l 

= C f E {L j [AL]L N -^}c (8) 

0<j<N-l 

= c f [AL N ~ l - L N ~ l A\ c = (9) 

By Jacobi's relation for commutators, |7 n ,7 m | is a constant of motion, and since this is a 
system that supports scattering, we see |7 n ,7 m J — > 0, and hence the system is completely 
integrable. 



III. THE TWO-BODY PROBLEM 



Having shown the system to be integrable, we then know the asymptotic wavefunction to 
be of the Bethe ansatz form, and the only input needed for the Bethe ansatz is the solution 
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to the 2-body problem. Such a solution is given in Landau and Lifshitz, Quantum Mechanics 
0. We summarize the results below. 

First, we discuss like particles. In terms of the relative coordinate r = x% — x%, the 
potential is s(s + l)/sinh 2 (r), and with the relative momentum k = (&2 — k\)/2, the wave 
function in the center of mass is given asymptotically as 

r s+1 r -> 0+, 

W(r) = < (10) 

[ e~ ikr + S{2k)e ikr r -» +00. 
The scattering amplitude S(k) is given by 

r(i + »fc/2)r(i + 3 -»fc/2) 

1 J T(l-ik/2)T(l+s + ik/2) { ' 

This scattering does not rearrange the particles. For bosons (fermions) the wavefunction 
must be (anti)symmetric, so the scattering amplitude for transmision will be ±S(k). In 
what follows, we will drop factors of —1 in the scattering amplitudes, assuming that they 
are taken care of by either the choice of statistics of the particles, the choice of quantum 
numbers as half-odd-integers, or the choice of number of particles as even or odd. 

Now, we discuss unlike particles, with the potential — s(s + l)/cosh 2 (r). The wave 
function in the center of mass is given asymptotically as 

e ikr + R(2k)e- ikr r -> -00, 
*(r)= (12) 

I T{2k)e lkr r -> +00. 

The reflection and transmission amplitudes are given as R(k) = S(k)r(k), T{k) = S(k)t(k), 
where 

/.s sin7rs 

sm7r(s + ik/2) 

smnik/2 iS 

t(k) = 7 14 

v ; sin7r(s + ifc/2) v 1 

Also present are bound states, labeled by an index m = 1,2,..., M, according to in- 
creasing energy. Thus the parity of the states is (— l) m_1 , and even parity states have spin 
0, while odd parity states have spin 1. Bound states appear as poles of the reflection and 



transmission amplitudes, R{k\ — k 2 ) and T{k\ — k 2 ) on the positive half of the imaginary axis, 
given by k\ >2 = k ± m, k > 0. The momentum and energy of such a bound state is P = 2k, 
and E = k 2 — k 2 . From the particular form of the reflection and transmission amplitudes, 
we find K, m = s + 1 — m, where m — 1,2, ... , M, and M(s) is the largest integer less than 
s. There are no bound states for < s < —1. Threshhold values of s are s = 0, 1, 2, . . ., 
and at these values, the reflection amplitude vanishes. At the bound state poles we also find 
r(2iK m )/t(2iK m ) = (— l) m_1 . We call the bound states pairs. 



IV. YANG-BAXTER EQUATIONS 

We know the Yang-Baxter equations must hold, and we can verify this explicitly. For 
a two-component system the Yang-Baxter equations are equivalent to r 2 = r 3 r\ + t 3 r 2 t\ 
and r 3 t 2 = r 3 ti + t 3 r 2 r 1 , where n = r(ki - k 2 ),r 2 = r(ki - k 3 ),r 3 = r(k 2 - fc 3 ) etc. for tj. 
A degenerate situation occurs at a pole in r 3 and t 3 , when k 2 — k 3 = 2m m . There, since 
r 3 /t 3 = (-l)" 1 ^ 1 , the equations become = r x + {-\) m - x r 2 t x and t 2 = t x + (-l) m - 1 r 2 r 1 , 
where r 2j i = r(k ± iK m ), and etc. for tj. These relationships will be important when we 
calculate phase shifts. 



V. PHASE SHIFTS 



If a particle of type m passes through a particle of type ml ', without reflection, then we 
have a scattering amplitude exp[— i0 mm '(ki — k 2 )} , and a corresponding phase shift 9 mm '(k). 
These phase shifts are at the heart of the Bethe ansatz. Let us label the unbound particle 
by m = 0. Then we have found 

~T(l+ik/2)T(l + s-ik/2) 



6 00 (k) = Hog 



(15) 



Y{l-ik/2)Y{\ + s + ik/2) 

As we explained, we will not include factors of —1 in the scattering amplitudes, so #oo(0) = 
In general 9 mm '(k) = —0 mm r(—k) = 9 m ' m (k), and we will find that we always have 9 mm '(0) = 
0. 
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Now, consider the scattering of a particle k\ on a pair of two particles with momenta 
k 2 ± iK m . Let k = ki — k 2 . Then using the degenerate Yang-Baxter equations, we find for 
the scattering amplitude 

S(k + iK rn )S(k - iK m )t(k + in m ) = exp[-i9 0m (k)]. (16) 

Using the explicit forms, we can verify that # 0m (/c) is real for k real. 

Finally, we view the scattering of a pair from a pair as the scattering of two particles 
with momenta ki±iK m from a pair with fc 2 ±m m /. This gives us a net phase shift mm >(k) = 
Qom'(k — iK m ) +9om'(k + iK, m ). Again, using the explicit forms, we can verify that mm >(k) is 
real for k real, and symmetric in m,m'. 

To summarize: We have iVf particles with a = +1, and N± with a — — 1, for a total 
of N = iVf + Ni. Let us assume iVf > A 7 "^. Further, pairs of up-down spins bind into a 
variety of bound states, or pairs, labeled by m, m — 1, . . . , M(s). Let there be N m of each 
type. Then the number of unbound particles is N — N — 2 Si< m <M We will call these 
simply particles from now on. They would correspond to spinons/ions in the spin/charge 
picture. Of these particles, we have N_i with spin down; let us call them spin waves. Clearly 
N-i =Nl- Ei< m <M N m , and N_, < N /2. 

We still must treat the dynamics of the spin waves, but since they are not "real" particles, 
but only correlations in the quantum numbers of particles, they have no momentum or energy 
directly. Thus, defining 

0, m — —1, 

1, m = 0, (17) 

2, m= 1,2, ...,M(s) 
then we can write the momentum and energy as 

P= E V m J2 k m, (18) 

-l<m<M km 

E = \ E VmJ2 k2 m- E N m^m, (19) 
-l<m<M k m l<m<M 

(20) 



7 



VI. SPIN WAVES 



Since particle-pair and pair-pair pass through one another with only a phase shift and 
no reflection, their interaction is in some sense trivial. However, particles do scatter from 
particles with reflection, and their interaction is not trivial. We write the asymptotic wave- 
function explicitly in the Bethe ansatz form, and for now consider only the N Q particles. We 
use the spin language, so u z {j) = ±1 according to whether the jth particle in the ordering 
X\ < ... < Xn has spin up or down. A choice for all cr z (j) we denote simply by a. Then 
asymptotically the wave function is given by 



#(x\a) -> A ( p \ a ) ex P 
p 



2J Xjkpj 

l<j<N 



(21) 



The summation is over all the N \ permutations of the momenta. We arrange the A(P\a) 
for fixed P as a column vector £(P). Then the Yang-Baxter equations ensure that we can 
find a consistent set of amplitudes A(P\a), by finding the simultaneous eigenvector of the 
N equations 



e 



IJ S(kj - k n )Xjj-i ■ ■ ■ X jA X j>No ■ ■ ■ X jd+1 £(I) = ai) (22) 

l<n<Af 

In this equation, the X^ n are operators given as 

l + t{k)^ l-t(k) ... 

+ T -^T [°x(j)vx{n) + (Ty(j)cTy(n)] , (23) 

where k = kp. — kp . 

These eigenvalue equations can in turn be solved by a Bethe ansatz for the N — 1 
overturned spins — the spin waves — on a lattice of N particles. These equations can be 
solved either: (i) directly, by the methods of Yang || ; (ii) with commuting transfer matrices, 
by the methods of Baxter H]; or (iii) by quantum inverse scattering methods of Faddeev and 



Takhtajan PH |. We are not aware that these equations have appeared before in the solution 
of a quantum many-body problem, although the low-density case has often appeared, for 
instance first in Yang's original solution for 5-function fermions. 
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The solution is sufficiently technical that we postpone discussion to a later publication. 
However, the result has many interesting physical consequences, and that is what we want 
to discuss in this letter. One finds for the eigenvalues of the previous equations, 



J n < No 3 ijjk^ + i(kj - X q )}/2 



(24) 



In this equation, the A's are the momenta of the spin waves, and are determined from the 
equation 

yr smir[s + i(\j- \ q )/2] yr sin tt [s -i(Xj - k n )}/2 = 1 

We now have our two final phase shifts, for particle-spin wave and spin wave-spin wave 
scattering: 



= Hot 



sm7r \s 



ik}/2 



sin7r[s — ik]/2 
sin7r[s + ik/2] 



(26) 
(27) 



sin7r[s — ik/2] 

As noted, there is no phase shift for spin wave-pair scattering. 

At the threshhold values s = integer, the spin wave modes uncouple completely from the 
particles, and thus from the system, since they contribute no energy or momentum directly. 
In this case, we have the very high degeneracies found in the 1/r 2 lattice systems [11,12 
and for the same reason — an absence of reflection. 



VII. THE SOLUTION 

Let us now impose periodic boundary conditions and take any particle, pair or spin wave 
around a ring of large circumference. Along the way, it suffers a phase change as it scatters 
from every other particle, pair or spin wave, plus a phase change of PL, where P = rjk is its 
own momentum. Periodicity requires that this phase change be an integer multiple of 2tt, 
the integer being the quantum number. We write this statement as coupled equations in a 
rather symbolic form: 
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LTj rn k m 27T_/ m (/c m ) -|- ^ ^ . @m,m' ( k m ^m') 



m = 



-1,0,1,. ..,M. 



(28) 



-l<m'<M fc' , 



Here the I m (k m ) are the quantum numbers, the only subtlety being that for the spin waves, 
J_i ranges only over 1, . . . , Nq. 

VIII. RESULTS FOR ZERO TEMPERATURE AND ZERO SPIN/CHARGE 

In this letter, we give explicit results for the ground state and low-lying states when 
iVf = iVj, which we call the spin/charge zero-sector. This certainly is the most interesting 
case, since all singularities in the (N^, iVj_) ground state phase diagram occur for iVf = iVj_. 
In fact, as we shall see, for s > 0, the chemical potential has a discontinuity across the 
line iVf = iV|, and thus the system is an antiferromagnet /insulator, although not of the 
Neel/Mott type. For — 1 < s < 0, there is a weak singularity at iVj = iVj_, without a 
discontinuity in the chemical potential. 

For s > 0, which we call the attractive case, the ground state consists of a spin fluid 
of type m — 1, and thus spin 0. This is the bound state with lowest binding energy, when 
K — s, and so P(k) = 2k and E(k) = k 2 — s 2 . In the ground state, the fc's for the pairs 
distribute themselves densely with a density p(/c), between limits ±B, normalized so that 




(29) 



The energy and momentum are given by 




(30) 



(31) 



The integral equation which determines p(k) is 




(32) 



The kernel of the equation, 9' n (k), is the derivative of the phase shift. In figure (]!]) we show 
Eq/L versus N/L for selected values of s = 1/2, 1, 3/2. 
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Having determined the ground state properties of the system, we now determine the 
low-energy excited states. They are given by the following: (i) Remove a pair from the 
ground state distribution, and place it outside the limits; we call this creating a pair-hole 
and a pair, and it gives a two parameter continuum, (ii) Break a pair, to give two particles, 
one spin up and the other spin down; this also gives a two parameter continuum, (iii) Excite 
a pair into a higher energy bound state, if allowed; these we call excitons, and they have 
single parameter dispersion relations. (Away from the zero-sector, we can have in addition 
spin waves. These will be important for s < 0.) 



By the techniques of Yang and Yang ||13|| , the dispersion relations are given parametrically 

by 



m 



B 

Vmkm ~ / @ml(km ~ k)p(k)dk 



B 



Tj m k. 



2 



1 f B 

7T / e 'mi( k m - k)e(k)dk 

Z7T J-B 



(33) 



(34) 



2 

Here e(k) is the solution to the integral equation 

k 2 - s 2 = in = e(Jfe) + [ B eUk - k')e(k')dk'. (35) 

2,71 J-B 

The chemical potential ii\ is the chemical potential for pairs, given by dE /dNi. The results 
are shown in figure (@) for s = 3/2, B = 3/2, d = N/L = 0.943, E /L = -0.691, m = 1.215. 
The gap for the creation of two particles is AE = 1.170, and is equal to the discontinuity of 
the chemical potential across the line — iVj_. The exciton with m = 2 is the only exciton 
allowed at this value of s, and has a gap of AE = 1.017. 

For > s > — 1, in the zero-sector, we have two coupled equations for iV particles and 
N/2 spin waves. However, in the zero-sector, the limits of the spin wave distribution are 
±oo. Thus we can solve by Fourier transforms for the spin wave distribution in terms of the 
particle distribution, and then substitute this into the particle equation, giving us a single 
integral equation for the distribution of particles p(k): 



1- = p(k) + ^( B Q\k - k')p{k')dk'. (36) 

Z7T Z7T J-B 
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Here the kernel 9' 



(k) is given as 



oo 



dte 



ikt 



sinh 3 t(l + s)coshts 



(37) 



— oo 



In figure (|1|) , we show E /L versus Nj L for s 



1/2. The excited states in the zero-sector 



are given by: (i) Remove a particle from the ground state distribution, and place it outside 
the limits; we call this creating a hole and a particle, and it gives a two parameter continuum, 
(ii) Remove a spin wave from the ground state distribution, and place it on the line with 
imaginary part equal to i; we call this creating two spin waves, one with spin up and the 
other with spin down. It gives a two parameter continuum, familiar from the Heisenberg- 
Ising model. The results are shown in figure (^|), for s = —1/2, B = 1, d = N/L = 0.600, 



Finally, we remark that all thermodynamics can be explicitly calculated, since there are 
no ambiguities with counting states, or difficulties with strings of length greater than two. 
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Eq/L = 0.094, fi = 0.374. 
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FIGURES 



FIG. 1. Ground state energy per unit length Eq/L versus density N/L for s = —1/2, 1/2, 1,3/2. 

FIG. 2. Energy above the ground state energy versus momentum (dispersion relations) for the 
low lying excitations when s = 3/2 and density N/L = 0.943. 

FIG. 3. Energy above the ground state energy versus momentum (dispersion relations) for the 
low lying excitations when s = —1/2 and density N/L = 0.600. 
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